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A simple non-equilibrium thermodynamics is developed and a particular example is studied. The 
theory is formulated to describe a viscoelastic fluid, capable of finite deformation, which need not be 
locally in or near a state of thermodynamic equilibrium. This fluid may support shear stresses only 
when away from local thermodynamic equilibrium. A notion of time-temperature superposition is 
contained in the formulation of the constitutive equations. Conservation of energy is obeyed and the 
second law of thermodynamics is satisfied as a consequence of simple requirements on the constitutive 
relations. In an adiabatic isochoric motion the temperature increases when work is done on the ma- 
terial and decreases when the material does work. For given volume and temperature, entropy de- 
creases when the material is deformed from equilibrium. It is shown in what general way viscosity 
depends upon temperature. For infinitesimal strain, the special form of the stress-strain relations 
are derived in order to determine how temperature and time-temperature superposition enter in this 
case. 



1. Introduction 

In 1963 we established that a fluid theory gives the 
correct stress relaxation response for a large variety 
of materials, some of which, such as vulcanized butyl 
rubber, are usually treated as solids [l]. 1 We started 
then to investigate whether or not it would be fruitful 
to regard all matter as fluid. Subsequent experimental 
work [2] has indicated not only that the fluid point of 
view can be productive in a practical sense, but that 
the incompressible elastic fluid theory of our 1963 
article is also very promising. It is thus a natural 
step, at this juncture, to endeavor to extend our theo- 
retical understanding by attempting to create a thermo- 
dynamic theory, while deliberately adhering to the 
fluid point of view and bearing in mind that the in- 
compressible elastic fluid should represent an iso- 
thermal limit. We believe that the fluid notion is the 
key to a simple thermodynamics just as it was the key 
to the simple isothermal mechanical theory. 

We cannot reconcile the notion of a fluid with the 
notion of equilibrium with shear stresses present. 
Consequently from this viewpoint, even materials 
usually considered to be solid, such as vulcanized 
rubber, cannot be considered in equilibrium while 
supporting shear. If we are to describe shear, we are 
therefore left with no alternative other than to write 
a non-equilibrium thermodynamics. 

We offer a caveat. Because the non-equilibrium 
concept is essential to make the fluid point of view 
hold together, the reader who reflects on solidlike 
behavior of our fluid must not forget that we have 
given up thinking of a sheared state in terms of 
classical equilibrium thermodynamics. Otherwise 
he will arrive at paradoxes. 



1 Figures in brackets indicate the literature references <»n page 1 12. 



The problem of how to generalize equilibrium 
thermodynamics to include non-equilibrium relations 
is resolved by adjoining to the usual state variables a 
quantity which has dimensions of specific entropy and 
which depends upon the past history of the material. 2 
This quantity is the seat of irreversibility. With it we 
generate a statement of the second law which is 
meaningful for non-equilibrium. With it we also 
arrive at a thermodynamic compressible fluid theory 
which allows solidlike behavior, and which is con- 
sistent with the incompressible elastic fluid theory 
cited above. 

The basic definitions, assumptions, and equations 
which constitute the thermodynamic theory of perfect 
elastic fluids are given in section 2. Central to the 
understanding of these is the idea of local thermo- 
dynamic equilibrium or non-equilibrium and the rela- 
tions which hold in each case. Section 3 is devoted 
to the manner in which the notion of time-temperature 
superposition enters the theory. It is suggested that 
the reader who is anxious to get on to the questions of 
energy conservation and dissipation may skip or skim 
through section 3 in his first reading. In section 4 we 
write down the form which the equation of conservation 
of energy takes in the context of a perfect elastic fluid. 
Our statement of fading memory, together with non- 
negative heat conductivity, is shown to imply con- 
sistency with the second law of thermodynamics in 
section 5. Also for given specific volume and tempera- 
ture, the entropy will be less in local non-equilibrium 
than in local equilibrium. In section 6 we show that 
positive heat capacity at constant deformation implies 
that temperature is rising when work is being done on 
the material and falling when work is being done by 
the material isochorically and adiabatically. Upper 
bounds are obtained on the amount of work available 
by allowing a material to recover either adiabatically 
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or isothermally. In section 7 we derive a relation 
between temperature and viscosity which has been 
observed experimentally. And in the final section, 
section 8, we obtain the stress-strain relations for in- 
finitesimal strain. 



2. The Caloric Equation of State 

In the theory of inviscid gases, the caloric equation 
of state relates specific internal energy € to specific 
volume v and specific entropy s through an expression 
of the form 



€=f(v, s), 



(2.1) 



A stipulation of the form of the function / suffices to 
characterize the gas. The functions relating local 
pressure p and local absolute temperature T to specific 
volume and specific entropy are given by 



p = -f v (v,s) = -£ 



T=f 8 (v, s) = 



ds ' 



(2.2) 



(2.3) 



so that one has the familiar relation 
de = Tds— pdv. 

Although a given volume of gas need not always be 
in equilibrium because quantities such as temperature 
or pressure vary from point to point, yet at each point 
local thermodynamic equilibrium is assumed always 
to hold in the sense that (2.1), (2.2), and (2.3) hold 
pointwise. When we turn to viscoelastic behavior, 
we take the point of view that such behavior cannot 
be explained by local thermodynamic equilibrium, 
and that a fluid in thermodynamic equilibrium cannot 
support shear stress. Thus we set ourselves the task 
of creating a theory which deals with fluids which are 
not in local thermodynamic equilibrium. We do 
this with caution, using equilibrium thermodynamic 
theory (theory of gases) as a basic guide, and departing 
conceptually from this theory in as little measure as 
we feel will allow us to explain observed viscoelastic 
phenomena. The theoretical material thus described 
will be called a perfect elastic fluid. 

As will be made clear below, we incorporate into 
our theory the following notions 

1° There exists locally a specific internal energy e 
and a specific entropy s. 

2° There exists locally a non-negative quantity 2 
having dimensions of entropy per unit mass, which 
vanishes at thermodynamic equilibrium. Indeed, 
local equilibrium at a point prevails if and only if 
2 = at that point. 

3° In equilibrium a caloric equation of state of the 
form (2.1) holds. Whether or not local equilibrium 



prevails, the caloric equation of state reads 



(2.4) 



4° There is a local absolute temperature T which 
satisfies 



T=Mv 9 s + I). 



(2.5) 



5° The stress is given by a hydrostatic pressure 
-pdij plus terms that vanish at equilibrium where 



-/>,* + 2). 



(2.6) 



In developing our set of notions into a theory, we 
have the task of relating 2 and stress to other quan- 
tities. 3 To this end, we now discuss some kinematics. 

We use Cartesian tensor notation, letting the point 
x be given by Cartesian coordinates x%, i = 1, 2, 3. For 
a given motion, the position of a particle X at time 
t, x(£), can be expressed in terms of its position at 
time r, x(t), by a function 



Xi(t)=Xi[lL(T),t,T], 



(2.7) 



from which we may define the relative deformation 
gradients 



Xik(t, t) = 



dxj(t) 
dx k (r) 



(2.8) 



By a configuration, we simply mean an assignment 
of position x\ to each particle X. 

The Cauchy-Green relative strain tensor Cij{t, t) is 
defined by 

Cij{t, t) = x ki (t 9 T)x kj (t, t), (2.9) 

and the Green-St. Venant relative strain tensor is 

Eij(t, ^^[dj^^-dij]. (2.10) 

For the matrix of Eki{t, r) and Cki(t, t) we write E(£, r) 
and C(£, r) respectively. 

We return to the question of specifying 2 and stress 
away from local thermodynamic equilibrium, and 
proceed by assuming the existence of a function 

S[E(f,T),t] 

which is defined for all t 3= and depends on a strain 
E(£, t) only through 

/i(f,T)=x?+xs+x| 

/ 2 (£,r) = A*At + \tXl+X§\! 
/ 3 (f , t) = X?X1X| 



2 In 1950, P. W. Bridgman published a very perceptive article [3] in which he speculated 
on some of the properties which a non-equilibrium theory of this type would have to exhibit. 
We have here a concrete example of such a theory which meets Bridgman's requirements. 



''The incompressible elastic fluid [lj provided a key as to how to write the expression 
for stress. 
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where Af =Af(f, t), i= 1, 2, 3, are the principal values 

afCtf.r). 

We write 



S«[E(f,T) f t]»^ 

as 



Sr[E(f , r), t] 



ar 



(2.11) 
(2.12) 



these notations applying no matter what name is given 
to the arguments of S. 

Furthermore, we shall assume the existence of a 
material property br > 0, which scales time intervals 
in a manner to be made more precise in section 3 
below. The quantity b T depends on absolute tempera- 
ture T and has the value unity at a given reference 
temperature T . We now assert that the value of 
2 at present time t depends upon the past history of 
deformation and temperature through 4 



2«=fs E(*,t),£m0# 



br(r)dT, (2.13) 



where bjir) means the value of br at time r. 
we specify that stress cry at time t satisfies 5 



Further 



Vi, 



pSij 



+ 



J -c 



Xik(t, r)xjAt, r)S k i 



E(t 



> * r 



biOW 



brir)dT, 



where p=v l is density, Ski is defined by (2.11) and 
P by (2.6). 

(In (2.14), as subsequently, we do not mention ex- 
plicitly that quantities such as cr/j, p, 2, and T are 
evaluated at present time t when the context makes 
this clear.) 

We note that s + 2 may be eliminated between (2.5) 
and (2.6) to express p as a function of specific volume 
and local temperature, say 



p=g(v, T), 
whence (2.14) becomes 

cnj = -g{v, T)Sij 



•/; 



+ pT Xik(t, T)Xji(t, T)S k 



E(t 



,r),/; 



Mf)<# 



(2.15) 

(2.16) 
bT(r)dr. 



The expression (2.14) bears a resemblance to the 
constitutive equation of an incompressible elastic 
fluid given by Bernstein, Kearsley, and Zapas, which 
is a corresponding isothermal incompressible theory. 



4 If in further discussions it should prove desirable to have a name for 2, we suggest 
calling it entaxy. 

5 In equation (3.6) of reference [1], the potential U of the isothermal theory of the in- 
compressible fluid replaces pTS in eq (2.14) of the present article. 



Indeed, for constant temperature 



/: 



b T (£)d{ = (t-T)b T , 



in which case the stress is given by (2.14) as if each 
configuration that the material assumed in the past 
contributed to the stress as the preferred configuration 
of an elastic material, each contribution being weighted 
according to the time elapsed, added up, and modified 
by the addition of an appropriate hydrostatic pressure. 
Since the role that S plays is analogous to that of an 
elastic strain energy, it is reasonable to assume that 



S[E(«, r), f]?0 



(2.17) 



for all £ 52 and all values of E(£, r) for which S is 
defined, and that 



S[0,f]=S[E(*,*),f]=0 



(2.18) 



for all £ ^ 0, where is the zero matrix. 

From (2.17), it is clear that whenever S vanishes it 
assumes a minimum value, whence 

S k i[E, f]=0 whenever S[E, f] = 0. (2.19) 



(2.14) In particular, (2.18) and (2.19) yield 



S*/[0,f]=0 



(2.20) 



for all f ^ 0. 

Our theory has been set up to be consistent with 
1°, 2°, 3°, and 4°. To show that it is also consistent 
with 5°, we have but to establish that the integral terms 
in (2.14) vanish when local thermodynamic equilib- 
rium prevails, i.e., when 2 = 0. But if 2 = 0, we see 
from (2.13), (2.17), and b T > 0, that 



s[e&t),£«0#] 



vanishes for all r^t. Thus, (2.19) assures us that 
Ski also vanishes for all r ^ t and hence, from (2.14), 
<Jij= — p8ij, which completes our argument that 5° 
holds. 

In addition to (2.17) and (2.18), we assume fading 
memory in the form 



S t *S 0. 



(2.21) 



For an interpretation of the motivation for this assump- 
tion, let us consider the case in which temperature, 
and hence br, is constant, so that 



/: 



67<f)df = U-T)6r 



is proportional to the time interval between t and r. 
Consider, then, a given value of relative strain E, and 
consider the contribution to the internal energy which 
arises if the relative strain computed from the con- 
figurations at present time t and past time r takes 
this given value. Then, in view of (2.4) and (2.13), 
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(2.21) says that this contribution cannot increase for 
larger values of the time interval t — j. 

If up to and including present time the material has 
been motionless, or even just undeformed, i.e., 
E(*, t) = for r ^ t, then (2.13) and (2.18) imply that 2 
vanishes and thus local thermodynamic equilibrium 
prevails. But, furthermore, we can show that if our 
material is kept at constant deformation and at con- 
stant temperature from some time U onward, it will tend 
toward local thermodynamic equilibrium, i.e. If 
E(t, t) = and T(t) = T(ti) for all t, r ^ ti, then 

limS(t) = 0, 

where S(t) decreases monotonic ally when t > ti. 

Appendix A is devoted to the proof of this statement. 
Example: For an ideal gas 



as a standard or laboratory clock, which is not tempera- 
ture sensitive. 

If at time t = 0, the temperature sensitive material 
clock reads 6 , then at time t it reads 6{t) given by 



0(t) = Mf)d?+0o' 



(3.1) 



At a constant temperature T, the material clock 
will advance 



e = D v (i-y) exp^-, 
Lv 

where D > 0, y > 1 and CV > are constants. Gen- 
eralizing this equation slightly to allow a different de- 
pendence of e on volume, we may consider 



e = Y(v) exp 77- 



(2.22) 



where Y(v) is some given function of v. Starting 
with (2.22) in place of (2.1), we may form an elastic 
fluid for which (2.4) becomes 



e = Y(v) exp — (s -f 2), 



(2.23) 



where 2 is given by an expression of the form (2.13). 
Then from (2.23) and (2.5) we obtain 



6 = CvT. 

Also (2.15) becomes 

p=-[\nY(v)]'C v T, 
whence (2.16) reads 



(r ij =[]nY(v)YCvT8i. 



■ p t[ 

J-cc 



Xik(t, r)xji(t, r)S k i 



E{t 



■< 



H&di 



(2.25) 



(2.26) 
br(r)dT 



A perfect elastic fluid for which (2.23) holds may be 
called an ideal elastic fluid. 

3. Time-Temperature Superposition 

The idea of time-temperature superposition which 
we introduce is based on the notion of a clock be- 
longing to each particle of the material. Our assump- 
tion is that the material reads its own clock, which 
is temperature sensitive and runs at a rate br as fast 



0(t)-8(T) = (t-T)b T 



(3.2) 



when the laboratory clock advances t — r. At refer- 
ence temperature T Q the material clock and the lab- 
oratory clock run at the same rate, so that b T =l and 
6(t)-0(T) = t-r. 

Looking at (2.15) now, we see that p does not depend 
explicitly on the past history of deformation and tem- 
perature, but only on the present value of temperature 
and specific volume. Using (2.14) we see then, that 
[crij + p$ij] I pT depends entirely on the histories of 
the temperature and the strain relative to the present 
configuration. The concept of time-temperature 
superposition commonly used in linear viscoelasticity 
[4] can be stated in terms of these ideas as follows: 

Suppose that x = f(£, X) is a motion occurring at ref- 
erence temperature T with corresponding specific 
volume v(t, X) and stress (Tij(t, X) (where X is a particle 
label). If this motion is replaced by x = f(6^, X) at 
constant temperature T, with corresponding specific 
volume v(t, X) = v(brt, X) and corresponding density 
p(t, \) = p(b T t, X), then the corresponding stress 
&ij(t, X) will satisfy 



frijjt, X) + p(g, X)5i, _ o-jj(b T t, X) + p(b T t, X)8jj 



p(t 9 X)T 



p(b T t, X)T 



(3.3) 



(2.24) where 



and 



p(t,X) = g[v(t,X),To] 
p(t,X)=g[v(t,X)J]. 



It is a simple matter to generalize this statement to 
apply to situations in which the temperature varies 
with time. We replace each time increment dr by 
br(r)dT and restate the * principle as follows: Let 
x= f(t, X) be a given motion occurring at reference 
temperature T and corresponding specific volume 
v(t) with stress cn/t, X). Suppose that x = f(0(t), X), 
with corresponding specific volume v(t, X) = i;(0(t), X) 
and corresponding density p(t, X) = p(#(t), X) is a 
motion occurring at (possibly varying) temperature 
T(t) with corresponding bj(t), so that d(t) is given by 
(3.1). Let <Ty(t, X) be the stress corresponding to the 
second motion. Then 



& {j (t, X)+p(t, X)d ij ^ a ij (6(t\X)^p(e(t), X)dij 
p(t, X)T p(d(t), X)T 



(3.4) 
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inhere 



and 



p(t,X)=g(v(t,X),To) 



p(t,X)=g(v(t,X),T). 

This is our principle of time-temperature superposition. 
It is, of course, an assumption subject to experimental 
verification. 

To phrase our principle in a less precise but more 
intuitive way, suppose an observer sees a motion at 
temperature T . Suppose, then, that at temperature 
T(t) the observer, reading only the temperature sensi- 
tive clock, appears to see the same motion again. 
Then he will appear to see the same history of 

((Tij+pdijypT. 

It is now our job to establish that a perfect elastic 
fluid is consistent with this principle. To this end, 
consider the motion x = f(£, X) at temperature To. 
Then 

^)+f^ = r irdXkAttT)tt _ T]dr (3S) 

p{t)lo J-oc 

where 

7Tij[Xkl(t, T), t — r]=X ik (t, T)Xjit, T)Skl[E(t, T), t — t] 

(3.6) 

and we have dropped explicit mention of the particle X. 
If for the motion x(£) = f(£, X) we write for (2.7) 

x(*) = g[x(r), t, r], 

we shall get for the motion k{t) = x{6(t))={[6(t), X] 

M6(t)) = g[MO(T)),0(t\6(T)] 

k(t) = g[k(T),d(t\6(T)] 

whence the relative deformation gradients for the two 
motions, namely Xki[t, r] and Xki[t, r] are related by 



&u(t) + p(t) 
p(t)T(t) 



J —oo 



(0(t\e(O),e(t)-d({)]b T ({)dt. 



From (3.1) 

Q(t)-0({) = j^b r (Qd£. 
Putting (3.7) and (3.10) into (3.9) yields 



(3.9) 



(3.10) 



era 



= -p8 ij + pT\ TTij x k i(t,i)A 



Au(t,&,\ bjiQdC 



b T (m 



which, in view of (3.6), is exactly what (2.14) gives for 
cnj. Thus the perfect elastic fluid is consistent with 
our principle of time-temperature superposition. 

4. The First Law of Thermodynamics 

To write down the first law of thermodynamics it 
is convenient to distinguish two kinds of derivatives. 
The material derivative of a quantity <p will be written 

as <b or -p— . If we consider (pasa function of time and 
^ Dt r 

particle label (Lagrangian coordinates), the material 

derivative of <p is defined by 

The particle velocity v is given, for example, by 
Vk=Xk* The spatial derivatives of a quantity <p will 
be written <p f k- If we consider <p as a function of time 
and position (Eulerian coordinates), the spatial deriva- 
tives of (p are defined by 

<P,k = — W,x). 
dXk 

Conservation of mass, in the form of the equation of 
continuity, is written 



Xkl{t, T)=Xkl{d(t), 6(t)]. 

Now (3.4) and (3.5) yield, 

O-^Q+pU)_O-Q-(g(0)+pW0) 



p(t)T(t) 



p(0(t))T o 

Bit) 



p + pvk,k = 

\ ' or, in terms of specific volume 

V = Wk,k 

whence 



-I 



TTij[x k ,(9(t),T), d(t)-T]dr. (3.8) 



pv = v k , k = Vi,j8ij. 



(4.1) 



A/ri . , , r -iir >-/onx We write the first law of thermodynamics, which 

Making the change of variable from r to f , in (3.8), S8e8 conserV ation of energy, in the form 

where 



r=0(f), dr = b T d^ 



yields 



pe = ajjVij — q u 
where q is the heat flux vector. 



(4.2) 
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Now material differentiation of (2.4) yields, with whence, using (5.1) 
the aid of (2.5) and (2.6) 



ps- 



Differentiation of (2.13) yields 
t = S[E(t,t) 9 0]b 1 (t)+f t S k iE k Mr)dT 

J -oo 



(4.3) 



T 



T 2 



or 



pi + (f) f 0, 

n which implies that the Clausius-Duhem inequality 

+ b T (t) S t b T {T)dT. (4.4) holds locally as well as globally [5]. Thus with 

•^ _0 ° (5.1), the assumption of fading memory tells us that 



Now from (2.9) and (2.10) 



E kl -- 



! 



dvj{t) dxj(t) dxi(t) dvjjt) 



= 2«J 



dx k (r) dxi{r) dx k (r) dxi{r)_ 
dxjjt) dxjt) j_ dxjjt) dxjjt) 



+ - 



dx k {r) dxi(r) dx k {r) dxi(r 



(4.5) 



a perfect elastic fluid satisfies the second law of 
thermodynamics in the Clausius-Duhem form. 

Note that we may solve (2.3) for s, obtaining, for 
equilibrium, an expression of s as a function of v 
and T, say 

5 = h(v, T). 



Similarly, whether or not local thermodynamic 
where, of course, Vij = dvi(t)/dxj(t) may be expressed equilibrium holds, we may solve (2.5) for s + 2, ob- 
in terms of t and x(j) alone. taining 

Noting that, by virtue of (2.18), the first term on the 



right hand side of (4.4) vanishes, we obtain from (4.3), 
making use of (4.1), (4.4), and (4.5), 



b T (r) dr 



p6 = -pvij8ij + pTs 

+ pTvu J **(*, r)Xji(t, T)S ki E(t, r), J b T (g)d4 

+ pTbAt)f t St E(t 9 T) 9 f'ih{Z)dg b T (r)dr. (4.6) 

Using (2.14), we get for (4.6) 

p€ = (TijVij + pTs + pTbj(t)\ SMr)dr. (4.7) 

J — X 

Insertion of (4.7) into the first law of thermodynamics 
(4.2) yields 



s + l = h(v,T) 



s = h{v,T)-X. 



(5.3) 



(5.4) 



'/'■ 

J — oc 



pTs + q Ui = - pTb T (t) S t b T (r)dr 



In the next section we shall show a very important 
consequence of this equation. 



5. The Second Law of Thermodynamics 

The law that heat flows against the temperature 
gradient reads 



Recalling 2°, we see immediately from (5.4) the 
following: For a given value of v and T, the local 
entropy s has its maximum value at local thermo- 
dynamic equilibrium. Indeed, for a given configura- 
tion and T, the entropy has its maximum value when 
local thermodynamic equilibrium prevails. 

Suppose that from some time t\ onward the material 
is held at constant deformation, and that the tempera- 
ture is also maintained constant from t\ onwards. 
Then, as can be seen from appendix A, the value of 
local entropy s will tend towards its equilibrium value 
h(v, T) in a monotonically increasing manner. Further, 
(4.8) (5.3) tells us that s + X remains constant during this 
process, any decrease in 2 being exactly equal to an 
increase in s. 

Let us compare the values of local entropy for two 
cases, one in which local equilibrium prevails (S = 0) 
and another in which it does not prevail (X > 0). Then 
if the values of specific volume and temperature are 
the same in both cases, we shall find that the entropy 
is less in the latter case than in the former. In 
particular, at given v and 7 1 , the presence of shear 
stresses indicates that the entropy is strictly less than 
it would be at equilibrium. 



qtT,i ^ 0. 



(5.1) 



Coupled with the fact that b T > 0, (2.21) and (4.8) 
imply 



Example: For an ideal elastic fluid (2.23), we get for 
(5.4) 



pTs + qu^O, 



(5.2) 



s = Cv In ^77- — X. 
Y(i>) 



(5.5) 
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6. Temperature and Energy of Stretching Moreover under these circumstances, as we have 

just argued, 7"^ 0, whence 
Combining (2.4) and (5.3) yields r>(TY\ 

e=f(v,h(v,T)) = a(v,T\ (6.1) Dt 



where a is appropriately denned. Thus the specific 
internal energy may be expressed as a function of 
specific volume and local temperature. This relation 
holds whether or not local thermodynamic equi- 
librium prevails. 
We define C^ by 



r _da 
CA ~dT 



so that from (6.1), (6.2), and (4.2) one has 

da - n • 
p — v + pLbl =cnjVi i j — quu 

ov 

which, in view of (4. 1), could also be written 



pCbT = 



da 
<r,--S, 



Vi.j-Qi. 



(6.2) 



(6.3) 



(6.4) 



Now suppose that heat is flowing but no deformation 
is occurring so that the symmetric part of vtj vanishes, 
i.e., Vi, j H- Vj,i = 0. Then (6.4) becomes 

pC^T=—qi,i. 

Thus Ca is the specific heat capacity at constant 
deformation and we feel quite justified in assuming 
C A >0. 

On the other hand, suppose that the material is 
adiabatically undergoing isochoric deformation 
(sometimes called a shearing motion). Then (6.3) 
becomes 

pCbT=(TijVi,j, 

whence we deduce that in an adiabatic isochoric 
motion the temperature is rising or falling depending 
respectively on whether work is being done on or by 
the material. 

During an adiabatic deformation the specific work 
done is equal to the change in specific internal energy 
as implied by (4.2). Thus an upper bound on the 
specific amount of work which can be taken out of 
the material adiabatically is given by e provided that 
zero is a lower bound on e. Furthermore, when the 
material is giving up work adiabatically and iso- 
chorically, so that q\ i = 0, v = 0, (TijVi,j<0, then by 
(4.3) and (5.2), 

which together with (4.2) yields 

TX ^ VCTijVi^j. 



From this and n ^ it follows that TL forms an 
upper bound on the amount of specific work that may 
be extracted from the material in recovering adiabati- 
cally and isochorically. 

The specific Helmholtz free energy, F, may be 
defined by 

F^e-Ts. (6.5) 



From (6.1), (5.4), and (6.5), we get 



where 



<J>(v, T) = a(v, T)-Th(v, T). 



(6.6) 



It is clear from (6.6) that for given temperature and 
specific volume, or alternatively for given tempera- 
ture and configuration, F takes its minimum when 
local thermodynamic equilibrium prevails. Also, 
from appendix A, we see that for constant configura- 
tion and temperature, F will tend monotonically to 
its equilibrium value. 

From (6.5), (4.2), and (5.2) 

pF = pe — pTs — psT = <JijVi, j — psT 

— (pTs + qt, i) ^ <TijV},j — psT. (6.7) 

Thus, if <f>(t>, T) has a lower bound for each T when 
p ^ 0, then the difference between F and this lower 
hound constitutes an upper bound on the specific 
amount of work which may be recovered isothermally 
from the material. In any case, (6.6) and (6.7) tell us 
that TL gives an upper bound on the specific amount 
of work that may be recovered isochorically and 
isothermally from the material. 

If the material is deformed from equilibrium and 
then allowed to come back to its initial temperature 
and volume while still in shear, (6.1) tells us that there 
is no net change in internal energy, and thus, from 
(4.2), the work done on the material has all been con- 
verted into heat, which has been given off to the sur- 
roundings. If work is now extracted from the material, 
and then the material is allowed to return to its original 
volume and temperature, we similarly see that the 
material has effectively removed heat from its sur- 
roundings and converted it into work. When stresses 
relax to zero isothermally with no motion of the mate- 
rial, no energy changes or conversions take place. 
Rather, the material loses its potential to convert 
heat into work while in contact with isothermal sur- 
roundings. The loss of this ability can leave a net 
increase in the amount of heat that has flowed into the 
surroundings, giving rise to viscous heat production. 
Thus, in this theory, elastic behavior and viscous heat 
production are intimately related in contrast to their 
frequent introduction as separate phenomena [6] . 
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Example: For an ideal elastic fluid, (2.23), we have 
already obtained the special form of (6.1) namely 



Thus for this fluid 



CvT. 



C^ — Cv 



(6.8) 



(6.9) and 



Xikit, t) 



1 (t-r)y 



1 




and we see that the heat capacity at constant defor- 
mation is constant. Indeed, it is a simple exercise 
to establish that the most general perfect elastic fluid 
having constant C\ is 

e = Y(v) exp -|- [s + 2] + i//(v) (6.10) 



E(*,r) = 



{t~r)y 

(t-r)y 















where Y > and i// are functions of v and Cv > is a . 

constant. For (6.10), we also have (6.9). For (2.23) Thus we have ' for a sim P le shear ( 7J )' 

we get from (2.24) and (4.2) 



pCvT=(TijVi,j — qi,i, 



Xik{t 9 T)Xji(t, T)Ski \E(t, r), Mf)d? 

= <Pij (t-ryy,j^)di 



(7.2) 
whence we see that the temperature will be rising or 
falling when adiabatic work is done on or by the mate- for ap p r0 p ri ately defined functions a». i, /= 1, 2, 3. 
rial respectively. And this will hold whether or not Thus from (7 2) and (2 14) under isothermal cond itions, 
volume changes are occurring. 



For (2.23), (6.6) becomes 
F = C V T 

d>(v, T) = C V T 



^Mh 25 



l-ln 



CyT 
Y(v). 



<Ti2 = pT I <pi2 [(* - r)y, (^- r)h T ]b T dT. (7.3) 

Making the change of variable £ = (£ — t)6t in (7.3) 
yields 

o-i2 = pr (fvzlary^ £]dt; = pTH(ary) (7. 
Jo 



4) 



and thus, in light of (2.25), <£(u, T) has a lower bound 
for each given T and all p ^ if and only if Y(v) is 



where wr— l/br and // is the function defined by (7.4). 



bounded away from zero for all values of v for which The viscosity 7) is defined by 
[In Y(u)]'^0. In particular, Y(v) = [a(v-Vi) 2 -l-fe]- 1 , 
where a and b and V\ are positive constants, is so 
bounded. On the other hand, for an ideal gas 



0"12 



Y(v) = Dv 1 -y, D>0, y>l, 

we have Y(u)— »0 as v-^ 00 , whence 4>(v, 7) has no 
lower bound for any T. 



7. Viscosity-Temperature Relations 

The motion 

xi(t) = xi(t) + (t - T)yx 2 (r) 



whence, from (7.4) 

H(ary) 
V = PT—^-' 

The zero shear viscosity 170 is defined by 



(7.5) 



170 = hm 7] 

y-+0 



whence, from (7.5) 



r)o = pTa T lim ^ ]7 =pTa T H f (0). 



x 2 (t) = x 2 (r) 

X3(t) = X3(r) 



a T y 



► ary 



(7.1) Thus, from (7.5) and (7.6) 



with y constant, represents a simple shear. The 
relative deformation matrix is 



7] _ Hjary) 



770 aryH'(0) 



(7.6) 



(7.7) 
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i.e., 17/170 is a function of a T y alone. Now from (7.6) Using (2.15), let 

JhL- K(p , T)=g v (v , 7>o, 

where 



ar 



(8.7) 



gv = 



dv 



pTH'(0) 

whence, eliminating ar in (7.7) yields 

V_ H(7ioylpTH t (0)) _ L f7 ]o y\ 

170 VoylpT \ pTJ ' 

1 hen, to hrst order in v 
where L is an appropriately defined function. , T , 

Equation (7.8) is in agreement with experimental gfo T)=p + ' (v - v ) =p + K(p , 7>u-- (8.8) 
data reported by F. Bueche [7]. u o 

Substituting (8.3), (8.4), (8.5), (8.6), and (8.8) into 
(2.16), and dropping all terms of order greater than 
one in mj gives 

(Tij =—po&u — *(po, T)8tj 

+ dyTj <p ( J WWf) MrWrtre(0 

4- 7 e (0 £ j/r ( J W#) Mt¥t 
- rSy J' <p ( [ M0df) rre(T)Mr)rfr 



8. Infinitesimal Strain 

Let us expand S in powers of 

Ji = trE(t, r)=E kk (t, t)= id \-3) 
J 2 = trE 2 (t, T)=E ki (t, T)Ei k (t, t) 

= 1(/i-3) 2 + (/ 1 -3)- K/ 2 -3) 
J, = det EU, t) = i [/, + 7,-/2-1]. 



(8.1) 



Because S[0, f]=Sw[0, J] = 0, we get, to second 
order in E k i(t, r) , 



j[e(*,t),£ 



fc<f)# 



(8.9) 



where tre = ekk- If we separate the stress and strain- 
1/n \ 1/P \ mto pressure parts, trar = <Jkk, tre, and shear parts, 

= 2 <P \J r br ^ )d %) ^1 + 2 ♦ (I 6 ^*j ^ (8 ' 2) ^U = o- U -i^jtra, e u = e u - J8^re, we get 



where <p and 1// are appropriate functions. Thus ^(t) = tro-(t)-\-po 

S kl = (ptrE8 kt + ilfE k i (8.3) 



to first order in E k u 

In appendix B, we show that for infinitesimal strain, 
with small displacement u{t) from an undeformed 
configuration, we may take 



whe 



xik(t, t) = dik + Ui, k (t) — Ui, k {r) 
Eij(t, T) = e/jU)-e/j(T), 

eij=Hu>i,j + Ui,j]. 



= tre(t) 



fhere 



*(po, D - P K' (f M0#) b^dr] 
T P K' ( P bA&dA tre(r)bAr)dr (8.10) 



(8.4) 



K(t) = j"[3<p(T) + llKT)]dT 



and 



Assume now that the undeformed configuration is in 
equilibrium at temperature T and ambient pressure 
p with specific volume v — f(po, T) and density p — 



.5) d-ij(t) = - Te tj J"' C ( | t ' M£Mf) 67<T)rfr 

+ r /I G ' ( // Wf )df ) e (T)6 7 <T)rfT, (8.11) 



pipo, i ;. > 


^lanuaiu iiiiiiiiLesiiiiai siiam aiu 


uyhife yieius 

where 






p(*)=po[l-e*jfc(*)] 

v(t) = v [l+ekk(t)]. 


(8.6) 
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G(*) = J%(f)#. 



In the first integral of (8.11) we make the change of 
variable 



at 



■H 



bte)<% 



to get 

crtj = TGiOMt) + T J' G' ( V b^dA e^b^dr. 

(8.12) 

In changing variables, we have used the assumption 
that briO i s bounded away from zero so that £(£, — oo) 
is infinite. Then (8.12) follows upon observing that 
G(°°) is zero. Similarly we get from (8.10) 

tt = [-k(po, T) + TK(0)]tre(t) 

+ T T K' ( P b T (£)dA tre(r)b T (r)dT. (8.13) 

In particular, for isothermal stress relaxation, in 
which eij{t) = for t < 0, eij(t) = const, for t > 0, (8.12) 
and (8.13) yield 

fry — TG(b r t)eij 

7r=[-K(po, T) + TK(b T t)]tre. 

Example: For an ideal elastic fluid (2.25) gives us 

po = -[ln Y(v )]'CvT (8.14) 

whence v is a function of po/CvT. Moreover, from 
(8.7) 



K(p ,T)=-[ln Y{vo)]"CvTv . 



(8.15) 



Thus we may solve (8.14) for t» and insert the result 
into (8.15) to obtain an expression for K(p , T) of the 
form 

K(p , T)=-fi (r^j CvT, 

where (3 is an appropriately defined function. 

From this we see that the infinitesimal stress-strain 
relations for an ideal elastic fluid are given when the 
three functions G, K, /3, and the constant Cv are speci- 
fied. The relation for the pressure component spe- 
cializes for an ideal elastic fluid from (8.13) to 



9. Discussion 

In our presentation of the theory of perfect elastic 
fluids, we have concentrated on the mathematical con- 
sistency of the formulation, and have tried to keep our 
theory as simple as possible. Without essential 
changes in the theory, one could set 2 equal to a sum 
of terms, each similar to the right hand side of (2.13), 
but each involving a different br, or one could let b T 
depend upon specific volume v as well as temperature. 
Numerous other modifications are also possible. We 
feel that further investigation is needed before decid- 
ing whether such modifications are necessary to make 
our theory describe real materials. 

At present it appears to us that our simple formu- 
lation can describe elastomers in what rheologists call 
the rubbery and flow regions. The applicability of the 
present formulation to the glassy region cannot be de- 
cided without further investigation. 

Finally we should emphasize that a principal fea- 
ture of this paper is its new approach to non-equilibrium 
thermodynamics. It is not a theory of perturbation 
about equilibrium. Nor need there exist an equilib- 
rium path between two given equilibrium states. We 
feel that the essence of this approach should be appli- 
cable to any thermodynamic process. 
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Appendix A 

We are to show here that ifE(t , r) = and T(t) = T(h) 
for all t, r^ ti, then lim 2(j) — where 2<(t) decreases 
monotonically. 

Since E(£, r) = for t > t\, monotonicity follows 
from (4.4), (2.18), and (2.21). We have yet to show 
that lim £ = 0. 

T K' ( ( br(i)d^ J tre(r)bT(r)dT , Let 8 be any positive number. For t 5= t\ we have 



tre(t) 



but the relation for the 
simplify from (8.12). 



shear components does not q <c £ : 



= j^S^E(t u r),j' 



bii£)d{ 



6t(t)Jt. (A.l) 



112 



Because of convergence of the improper integral (A.l) 
for t = t\, there exists a t-i < t\ such that 



j' 2 S E(t uT ), j' H€)d4 



6r(T)dT<2 



(A.2) 



and because of (2.21) this must also hold for t^ t\. 

Digress now to consider a motion involving only a 
single step deformation at £ = 0, such that T=T . 
For such a motion £ becomes at zero time just after 
the step 



Appendix B 

We discuss here in a non-rigorous manner the re- 
duction of finite strain quantities to their infinitesimal 
strain counterparts. We suppose that a particle is 
displaced from its initial position Xi at time t by 
Ui(t , X). Thus the position of the particle X at time 
t, namely x\(t) is given by 



/: 



S[E(0, r\t-r]dT 



whence 
Thus 



where E = E(0, r) is constant. In order for this in- 
tegral to converge with St ^ 0, we must have 

limS[E,£|=0 

for each value of E, and also this convergence must be 
monotonic. Returning now to (A.l) we thus see that 
for each fixed T, t% ^ T ^ tu 



Xi(t)=Xi + iH(t,X) 
Xi(i) = Xi{r) + m(t, X) — ih(t, X). 

dX t 



(B.l) 



Xik(t, t)= S ik + 



dxkir) 



dw(t, X ) dUiJT, X ) 
dXi dXi 



E(t 



,t), Pmo*] 



tends monotonically to zero as t->°°, whence it must 
do so uniformly in t, t->^r^ t\. (Dini\s theorem.) 
Hence for large enough t 



V l Smt u r\ P bM)d^ 



br{T)dr<~ 



(A.3) 



Adding (A. 2) and (A.3) gives, with the use of (A.l), 
that for large enough t, we have |£| <8 whence our 
result is established. 
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= 8 jfc + /f\y [Ui,m(t) - »i,m(T)] (B.2) 
dxk(v dXi 



where m, m (t) = dw(t, X)ldx m (t) and m, m (r) = dui(r, X/ 
dx m (t). 

From (B.l) 



dX[ _ dui{t,X) 
— oik — -— 

BXk(T) 

dUmit, X) 



dxk(r) 
dx m (t) 



dXi 



— S w /-h 



dX, 



(B.3) 



(B.4) 



Putting (B.3) and (B.4) into (B.2) and retaining only 
terms up to second order in the derivatives of dis- 
placement yields 



Xik(t, t) = ; 8** 4- m,k(t) — UjjAr) 



(B.5) 



to first order. 

Inserting (B.5) into (2.9-10) and retaining only 
first order terms in the derivatives of displacement 
yields 

Eij{t,T) = eij{t)-eij{T) 

to first order, where ey is defined by (8.5). 
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